Energy conditions in modified Gauss-Bonnet gravity 
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In considering alternative higher-order gravity theories, one is liable to be motivated in pursuing 
models consistent and inspired by several candidates of a fundamental theory of quantum gravity. 
Indeed, motivations from string/M- theory predict that scalar field couplings with the Gauss- Bonnet 
invariant, G, are important in the appearance of non-singular early time cosmologies. In this work, 
we discuss the viability of an interesting alternative gravitational theory, namely, modified Gauss- 
Bonnet gravity or f(G) gravity. We consider specific realistic forms of /(G) analyzed in the literature 
that account for the late-time cosmic acceleration and that have been found to cure the finite-time 
future singularities present in the dark energy models. We present the general inequalities imposed 
by the energy conditions and use the recent estimated values of the Hubble, deceleration, jerk and 
snap parameters to examine the viability of the above-mentioned forms of f(G) imposed by the 
weak energy condition. 
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I. INTRODUCTION 

An interesting approach in explaining the late-time ac- 
celerated expansion of the Universe [lj is the possibility 
that at large scales Einstein's theory of General Relativ- 
ity may break down. In this context, the Einstein field 
equation was first derived from an action principle by 
Hilbert, by adopting a linear function of the scalar cur- 
vature in the gravitational Lagrangian density. However, 
one may generalize the latter approach by considering 
higher order curvature invariants in the gravitational La- 
grangian density @, E|. The motivation for this proce- 
dure consists in the analysis of strong gravitational fields 
near curvature singularities and in considering consistent 
candidates of a fundamental theory of quantum gravity. 
Indeed, string/M-theory predicts that scalar field cou- 
plings with the Gauss-Bonnet invariant G are important 
in the appearance of non-singular early time cosmologies. 
These motivations may also be considered in the context 
of the late-time cosmic acceleration 0-0] )■ 

An interesting alternative theory is modified Gauss- 
Bonnet gravity, or f(G) gravity, wheref (G) is a general 
function of the Gauss-Bonnet term [8l4ll| . Note that the 
linear Gauss-Bonnet term is a topological invariant and 



the variation of the density ^—gG in the action leads to 
a total divergence, and therefore does not contribute to 
the field equations. Specific realistic models of /(G) grav- 
ity were constructed to account for the late-time cosmic 
acceleration [III and it is these forms of f(G) that 
we consider in this work. The respective constraints of 
the parameters of the models were also analyzed in [l2j]. 
More specifically, in [l2j the four types of finite-time fu- 
ture singularities emerging in the late-time accelerating 
era were studied from /(G) gravity. It was shown that 
by taking into account higher-order curvature corrections 
the finite-time future singularities in /(G) gravity are 
cured. Therefore, it turns out that adding such a non- 
singular modified gravity to singular dark energy models 
makes the combined theory to be non-singular one as 
well. In this context, we further consider the constraints 
imposed by the energy conditions and verify whether the 



parameter range of the specific models considered in 12 1 
are consistent with the energy conditions. More specif- 
ically, we define generalized energy conditions for /(G) 
modified theories of gravity, and consider their realization 
for flat Friedmann cosmological models. In particular, we 
analyze whether the weak energy condition is satisfied by 
particular choices of /(G) which were advocated in Refs. 
[Til [l2| as leading to viable models. 
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The energy conditions are fundamental to the singu- 
larity theorems and theorems of classical black hole ther- 
modynamics (we refer the reader to [13j for more de- 
tails) . Note that the energy conditions are obtained when 
one refers back to the Raychaudhuri equation for expan- 
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sion, where the attractive character of gravity is reflected 
through the positivity condition, i.e., Ruuk^k" > 0, with 
R pv the Ricci tensor and fc M any null vector. Now, in 
general relativity, through the Einstein field equation one 
ends up with T^k^k" > 0, which is the null energy con- 
dition. In particular, the weak energy condition (WEC) 
assumes that the local energy density is positive and 
states that T^U^U" > 0, for all timelike vectors If 2 , 
where T^ v is the stress energy tensor (for a perfect fluid 
we have p > and p + p > 0). By continuity, the WEC 
implies the null energy condition (NEC), T^k^k" > 0, 
where fc M is a null vector (l3| . The energy conditions 
have been extensively analyzed in the literature, such as 
in the cosmology settings and f(R) gravity and we refer 
the reader to Refs. (l4lTl5j for more details. 

This paper is outlined in the following manner: In Sec- 
tion UH we present the gravitational field equations for 
modified Gauss Bonnet gravity, and in Section Mil we 
outline the respective inequalities from the energy condi- 
tions. In Section [TV] we consider specific forms of /(G), 
and analyze the constraints arising from the energy con- 
ditions. Finally, in Section[V]we present our conclusions. 
Throughout this work, we consider the following units 
c = Q = 1 (here Q is the Newtonian gravitational con- 
stant to distinguish it from the Gauss-Bonnet term, G). 



II. FIELD EQUATIONS OF /(G) MODIFIED 
GRAVITY 

An interesting alternative gravitational theory is modi- 
fied Gauss-Bonnet gravity, which is given by the following 
action: 

S = ^2 J d A x^[R + /(G)] + S M (g^ ,V) , (1) 

where the Gauss-Bonnet invariant is defined as 

G=R 2 - AR^R^ + R^R^ . (2) 

It is also important to note that in the matter action, 
matter is minimally coupled to the metric and not to 
the scalar field, making Gauss-Bonnet gravity a met- 
ric theory. Thus, using the diffeomorphism invar iance 
of Sm{9 ,JjV ,ip) yields the covariant conservation of the 
stress-energy tensor, V^Tf^ = 0. Modified /(G) grav- 
ity has been extensively analyzed in the literature and 
instead of reviewing all of its intricate details here, we 
refer the reader to p-Q; 

Now varying the action ([TJ with respect to the metric 
provides the following gravitational field equation 



~ 2 Rg ^ = K2T ^ at) + 5 + (~2RRw + ^ P Ru P - 2R/° T R vpaT + Ag a Pg^ R pavf} R p<T ) f {G) 
+2 [V^V„/'(G)] R - 2 9fW [□/'(<?)] R + 4 p/'(G)] - 4 [V p V p f(G)] R/ 
-4 [V p V u f(G)\ R/ + Ag pv [V P V CT f'(G)] R»° - 4 [V P V CT /'(G)] g ap g^R pa ^ , (3) 



where the prime denotes differentiation with respect to 
G. Note that is the covariant derivative opera- 
tor associated with g pu , □ = ^"V^V^ is the covari- 
ant dAlembertian, and T^ at ' is the contribution to the 
stress energy tensor from ordinary matter. 

In this paper, we consider the flat FRW space-time 
described by the metric 

ds 2 = -dt 2 +a 2 (i)dx 2 , (4) 

where a(t) is the scale factor. 

In the FRW background, and taking into account a 
perfect fluid equation of state for ordinary matter, it fol- 
lows that the field equations for /(G) gravity are given 
by 

2AH 3 f'{G) + 6H 2 + /(G) - G/'(G) = 2n 2 p , (5) 
8H 2 f'(G) + l6Hf'(G) (H + tf 2 ) + (AH + 6# 2 ) 
+/(G)-G/'(G) = -2 K 2 p, (6) 



where p and p are the energy density and pressure, re- 
spectively, and the overdot denotes a derivative with re- 
spect to the time coordinate, t. 
Moreover, we have 

R = 6 (2H 2 + Hj , (7) 
G = 2AH 2 (h 2 + ff) . (8) 

In the FRW background, the gravitational field equa- 
tions may be rewritten to take the following form 

p cS = \h 2 , PcS = -1^(2H + 3H 2 ) , (9) 

where p e s and p e g are the effective energy density and 
pressure, respectively, defined as 
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PcS 



2k 2 



-/(G) + 24ff 2 (h 2 + H\ /'(G) - 24 2 iJ 4 ^2F 2 + HH + 4H 2 H^j /"(G) 



+8H 4 H + H 2 'ii f"{G) + 8{2A) 2 H i (2B ~ 

where Eqs. JT])-© were used. 

We also present the following useful relationship 



2H Z + HH + AH 2 H 



QH A + 8HHH + 2AH Z H Z + 6H H 

2 



/"'(<?)} 



PcS + PcS 



P + P + 



96H 2 



K 

+24H 2 (2H 2 + HH + AH 2 ^ /'"(G) 



(qH z + 8HHH - 18H 2 H 2 + 3H 3 H - 4H 4 H + H 2 H^Jf"(G) 



as it will be used throughout the text in the context of the energy conditions. 



r 



(10) 



(11) 
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III. ENERGY CONDITIONS 



The energy conditions arise when one refers to the Ray- 
chaudhuri equation for the expansion, given by 



d6 



j pu a pu + u^" 



■Rio/ k, 



(13) 



where i? M „ is the Ricci tensor, and 9 , a pu and oj^ are, 
respectively, the expansion, shear and rotation associ- 
ated to the congruence defined by the null vector field 
k p . Note that the Raychaudhuri equation is a purely ge- 
ometric statement, and as such it makes no reference to 
any gravitational field equations. 

The shear is a "spatial" tensor with a 2 = u^a^ 1 ' > 0, 
thus from Raychaudhury's equation it is clear that for 
any hypersurface orthogonal congruences, which imposes 
(jjnv = 0, the condition for attractive gravity reduces to 
R^ u k^k v > 0. The latter inequality ensures that geodesic 



congruences focus within a finite value of the parameter 
labeling points on the geodesies. However, in general rel- 
ativity, through the Einstein field equation one can write 
the above condition in terms of the stress-energy tensor 
given by T ilv k p, k v > 0. In any other theory of gravity, 
one would require to know how one can replace us- 
ing the corresponding field equations. In particular, in a 
theory where we still have an Einstein-Hilbert term, the 
task of evaluating R^k^k" is trivial. However, in /(G) 
modified theories of gravity under consideration, things 
are not so straightforward. 

For convenience Eq. ([3]) may be written as the follow- 
ing effective gravitational field equation 
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where the effective stress-energy tensor is given by 



(14) 



T£ = k 2 T^ + -g^f(G) + {-2RR^ + AR W R V P - 2R/° T R upar + Ag ap g^R^ vl3 R pa )f'{G) 
+2 [V M V„/'(G)] R - 2 9fiu [□/' (G)] R + 4 [□/'((?)] - 4 [V P V^/'(G)] RJ> 

-4 \y P Vvf'(G)] R/ + [V p V ff /'(G)] R pa - 4 [V p V ff /'(G)] g ap g pa R^ . (15) 



In this context, the positivity condition, R^ u k p k u > 0, 
in the Raychaudhuri equation provides the following form 
for the null energy condition T°^k fJl k u > 0, through the 
modified gravitational field equation (fj"4|) . We also im- 
pose the condition T^h^k" > for ordinary matter. 
This is useful as applying local Lorentz transformations 
it is possible to show that the above condition implies 
that the energy density is positive in all local frames of 
reference. 



Taking into account that the Raychaudhuri equation 
holds for any geometrical theory of gravitation, we will 
maintain its physical motivation, namely, the focussing of 
geodesic congruences, along with the attractive charac- 
ter of the gravitational interaction to deduce the energy 
conditions in the context of /(G) modified gravity. To 
this end, using the modified (effective) gravitational field 
equations the energy conditions in this context are given 
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by 



NEC^Peff+Peff >0, 



(16) 



WEC p cS > and p cS + p cS > , (17) 



SEC 



Pes + 3f> e fi > and p eS + p cS > , (18) 



DEC 



Pes > and p e s ± p e s > , 



(19) 



where the notation NEC, WEC, SEC and DEC stand for 
the null, weak, strong and dominant energy conditions, 
respectively. 

Now, in standard mechanics terminology the first four 
time derivatives of position are referred to as velocity, 
acceleration, jerk and snap. In a cosmological setting, in 



addition to the Hubble parameter H = a/a, it is appro- 
priate to define the deceleration, jerk, and snap parame- 
ters as 

la l a l a' 

<l = - — -i J = TTr— . and s = T7Z— '( 20 ) 



H 2 



H 3 a ' 



I/ 4 



respectively. 

In terms of the these parameters, we consider the fol- 
lowing definitions 



H = -H 2 (l + q), 
H = H 3 (j+3q + 2) , 
H = H 4 (s-2j-5q-3) 



(21) 
(22) 
(23) 



respectively. 

Using the above definitions, then the energy conditions 
(fT5|) - (p^|) take the following respective forms 



96 

NEC : peff + Pes = P + P + p { -(6a 3 + 27a 2 + 21a + 8qj + 9j - s)f"(G)- 



24[4(a 2 + 2q + 1)H 2 + 2q 2 +7q + j + 4]/"'(G)} H 8 > , 



(24) 



WEC: p cS = p+—[-f(G)-24H 4 qf(G)-(24) 2 H s (2q 2 + 3q + j)f'(G)]>0, Pc s + Pc s>0, (25) 



SEC : p cff + 3 PcS = p + 3p+ [/(G) + 2AH A qf'(G) + 288 J ff 8 (-6g 3 - 23a 2 - 15a - 8qj - 7j + s)f"{G) 



+(24)(288)i/ i2 (2g 2 + 3g + j) 2 /"'(G)] > , p cS + PcS > , 



(26) 



DEC 



1 



PcS - PcS 



= P-P+ 72 [-/(G) - 24tf 4 g/'(G) - 96ff 8 (-6g 3 - 15g 2 - 3g - 8gj - 3j + s)/"(G) 



fc 2 

-(24)(96)ff 12 (2g 2 + 3g + i) 2 /"'(G)] > , 

I 



PeS + PeS > 0, p e fi > . 



(27) 



IV. VIABLE /(G) THEORIES USING THE 
ENERGY CONDITIONS 



Viable /(G) modified theories of gravity were used 
in [ll[ to account for the late-time cosmic acceleration. 
These latter models were studied in the context of curing 
the four types of finite-time future singularities emerging 
in the late-time accelerating era [12| . Indeed it was shown 
that by taking into account higher-order curvature cor- 
rections, in the context of /(G) gravity, the finite-time 
future singularities are cured. In this context, we fur- 
ther consider the constraints imposed by the energy con- 
ditions and verify whether the parameter range of the 
specific models considered in [12J are consistent with the 
energy conditions for flat Friedman cosmological models. 

Thus, we consider some specific forms of /(G), consid- 



ered in [n], [l2| given by 



h(G) 



ai G n + 61 



a 2 G" + 6 2 
/2(G) = a 3 G"(l + o 3 G m ). 



(28) 
(29) 



where a±, a 2 , 61, o 2 , 03, 63, n, and m are constants. In 
the following, we always assume n > 0. 

Note that the Gauss-Bonnet invariant, defined in Eq. 
(O, can be expressed as 



G = -24i? 4 g , 



(30) 



in terms of the Hubble and the deceleration parameters, 
respectively. 

As the inequalities imposed by the energy conditions 
in /(G) gravity are extremely lengthly, in the follow- 
ing analysis we only consider the WEC in exemplifying 
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the application of the energy conditions. We consider 
the following present-day values for the deceleration, 
jerk and snap parameters [1(1 E3: <7o = —0.81 ± 0.14, 
j = 2.16±g;?l, and s = -0.22+°^. 



uum, i.e., p = p = 0. The WEC constraints, i.e., p e ff > 
and p e ff + Pes > 0, are respectively given by 



A. Specific case: ,/i(G) = fig^fc 



In first place, we consider the specific case of Eq. (1281) . 
For simplicity in the examples analyzed we consider vac- 



-[ai(-24gi/ 4 )" + 6i][a 2 (-24giJ 4 )™ + b 2 ] + n(-24< Z i/ 4 )"(ai6 2 - a 2 b 1 ) + {24y?iH s (-2AqH*) n - 2 [a 2 (n + 1) x 



n(-24 9 i/ 4 )"(-6 1 a 2 + ai6 2 ){(6g 3 + 27q 2 + 21q + 8qj + 9j - s)[n{al{-2AqH A ) 2n - b 2 2 ) + ( a2 (-24 9 if 4 )" + 6 2 ) 2 ] 

+ {AH 2 + 8H 2 q + AH 2 q 2 + 2q 2 + 7q + j + 4)[4a 2 6 2 (l - n 2 )(-24giJ 4 )™ + a\(n 2 + 3n + 2){-2AqH 4 ) 2n 

+bl(n 2 -3n + 2)]q- 1 H- 4 }>0. (32) 




FIG. 1: The plots depict the weak energy condition for the specific form of /i(G) = '^ n ^' . The left plot corresponds to 
Pcff > 0; the right plot corresponds to p c fr + p c « > 0. We have considered the values a\ = —1, = —1, and 02 = 2. The plots 
show that the weak energy condition are satisfied for the parameter range considered. See the text for details. 



The constraints provided by the inequalities (|3ip - (|3"2")l 
are too complicated to find exact analytical expressions 
for the parameter ranges of the constants a\, a 2 , b\, & 2 , 
and n, so we consider specific values for some of the pa- 
rameters. In particular, we impose the following values 
di = —1, bi = —1, and a 2 = 2, and plot the WEC as a 
function of 6 2 and n, which is depicted in Fig. [TJ The 
latter does indeed prove that the specific form of fi(G) 
given by Eq. ([28]) considered in [12| is consistent with 
the WEC inequalities. 



B. Specific case: f 2 {G) = a 3 G n (l + b 3 G m ) 



We consider the specific realistic case of Eq. (|29]l an- 
alyzed in il2j which accounts for the late-time cosmic 
acceleration, and that cured the four types of finite-time 
future singularities emerging in the late-time accelerat- 
ing era, given by the following specific conditions n > 0, 
to < and n 1 and for several parameter ranges. 
Rather than exhaustively analyze all of the cases, we 
consider a specific case that does indeed prove that in 
addition to curing the finite-time future singularities is 
satisfies the weak energy condition. The latter parame- 



-0.00052 
-0.0005 



-0.0005 



FIG. 2: The plots depict the weak energy condition for the specific form of /2(G) = a?,G n {l + &3G" 1 ). The left plot corresponds 
to p e ff > 0; the right plot corresponds to p c g + p e ff > 0. The parameter range for this specific case is given: n > |,n 7^ l,m < 
0, n + m > 1, 0363 > 0. We have considered the specific values (n = 2.5, m = —1). See the text for details. 




FIG. 3: The plots depict the weak energy condition for the specific form of /2(G) = a3G n (l + 63G" 1 ). The left plot corresponds 
to p e ff > 0; the right plot corresponds to p c s + p c ff > 0. The parameter range for this specific case corresponds to: n > i,n =^ 
1, m < 0, 1 < n + m < 1, 0363 < 0. We have considered the specific values (n = 1.8, m = —1). See the text for details. 



ter range is given by following For the form of /2(G) considered by Eq. ([29]) . the 

WEC constraints, i.e., p e g > and p e g + p c s > 0, are 
n > 1/2, n + m > 1 and 0363 > , (33) given by 

n > 1/2, 2/3 < n + m < 1 and a 3 b 3 < , (34) 

n > 1/2, and n + m < 2/3 . (35) 



-a 3 {(-24 9 i/ 4 ) ,l [l + 6 3 (-24 9 i? 4 ) m ] + 2AH i {-2AqH i ) n - 1 [n + (n + m)b 3 (-24qH 4 ) m ] 
+24 2 H 8 {-24qH A ) n - 1 [n+ (n + m)6 3 (-24giJ 4 ) m ](2g 2 + 3q + j)} > 0, (36) 

a 3 (-24qH 4 ) n {(6q 3 + 27q 2 + 21q + 8qj - s)[n 2 - n + b 3 (-24qH 4 ) m {n 2 -n+ 2nm + m 2 - to)] 
+ (4H 2 + 8H 2 q + 4i/ 2 g 2 + 2q 2 + 7q + j + 4)[n 3 - 3n 2 + 2n + 6 3 (-24 9j ff 4 ) m (n 3 - 3n 2 + 3n 2 TO 
+2n - 6nm + Znm 2 + to 3 - 3m 2 + 2TO)]g" 1 J ff" 4 } > , (37) 

I 



respectively. 



As in the previous example, the constraints provided 
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FIG. 4: The plots depict the weak energy condition for the specific form of /2(G) = as,G n {l + 63G" 1 ). The left plot corresponds 
to p e fi > 0; the right plot corresponds to p c « + p c ff > 0. The parameter range for this specific case corresponds to: n > ^ 
1, m < 0, n + m < |. We have considered the specific values of n = 1.3, m = — 1. See the text for details. 



by the inequalities (|36t - (|37| are too complicated to find 
exact analytical expressions for the respective parameter 
ranges of the constants 0,3, 63, m, and n, so we consider 
specific values for the parameters. The parameter con- 
straint given Eq. (1331) . with the following specific values 
(n = 2.5, m = — 1) are depicted in Fig. [5] the constraints 
provided by Eq. (TM)) are depicted in Fig. [3] for the val- 
ues (n = 1.8, m = —1); and finally the constraints pre- 
sented by Eq. (|55l) are depicted in Fig. 2] for the values 
(n = 1.3, in = —1). The respective WEC conditions are 
then provided as a function of the parameters 03 and 63. 
As in the previous case, the plots depicted in Figs. [2j4] 
do indeed prove that the specific form of /2(G) given by 
Eq. (|2T)|) considered in [l2[ is consistent with the WEC. 



V. DISCUSSION AND FINAL REMARKS 

The standard model of cosmology is remarkably suc- 
cessful in accounting for the observed features of the Uni- 
verse. However, there remain a number of fundamental 
open questions at the foundation of the standard model. 
In particular, we lack a fundamental understanding of the 
acceleration of the late universe. Recent observations of 
supernovae, together with the WMAP and SDSS data, 
lead to the remarkable conclusion that our universe is 
not just expanding, but has begun to accelerate. One is 
liable to ask: What is the so-called 'dark energy' that is 
driving the acceleration of the universe? Is it a vacuum 
energy or a dynamical field ("quintessence")? Or is the 
acceleration due to infra-red modifications of Einstein's 
theory of General Relativity? How is structure forma- 
tion affected in these alternative scenarios? What will 
the outcome be of this acceleration for the future fate of 
the universe? 

The aspects of these fundamental questions whose res- 
olution is so important for theoretical cosmology, need 
to look beyond the standard theory of gravity. A very 
promising way to explain these major problems is to as- 



sume that at large scales Einstein's theory of General Rel- 
ativity breaks down, and a more general action describes 
the gravitational field. It is clear that these open ques- 
tions involve not only gravity, but also particle physics. 
String theory provides a synthesis of these two parts of 
physics and is widely believed to be moving towards a 
viable quantum gravity theory. Thus, in considering al- 
ternative higher-order gravity theories, one is liable to be 
motivated in pursuing models consistent and inspired by 
several candidates of a fundamental theory of quantum 
gravity. In this context, predictions of string/M-theory in 
the context of gravity-matter couplings, show that cou- 
plings of the scalar field with higher order curvature in- 
variants are important. In particular, a coupling of the 
scalar field with the Gauss-Bonnet invariant G are fun- 
damental in the appearance of non-singular early time 
cosmologies. 

In this work, we discussed the viability of an inter- 
esting alternative gravitational theory, namely, modified 
Gauss-Bonnet gravity or f(G) gravity. We considered 
specific realistic forms of f(G) analyzed in the literature 
that account for the late-time cosmic acceleration and 
that cured the finite-time future singularities [IJ (HI- 
The general inequalities imposed by the energy condi- 
tions were outlined and using the recent estimated values 
of the Hubble, deceleration, jerk and snap parameters we 
have shown the viability of the above-mentioned forms 
of f(G) imposed by the weak energy condition. More 
specifically, for simplicity we only examined the vacuum 
case for which p = p = 0, although this is not a phys- 
ically interesting the universe contains matter 
and radiation. However, this is easily corrected, since 
one can always add a positive energy density or pressure 
from matter and/or radiation to any model satisfying the 
WEC, and it will still satisfy the WEC. Thus, f(G) grav- 
ity with matter will also satisfy the WEC if the vacuum 
model does. 

However, as argued in [X 51 ] it is important to emphasize 
that although the energy conditions in modified theories 
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of gravity have a well-founded physical motivation, i.e., 
the attractive nature of gravity as outlined in Raychaud- 
huri's equation, the issue as to whether they should be 
applied to modified theories of gravity is an open ques- 
tion, which is ultimately related to the confrontation be- 
tween theory and observations. 



RGC grant of the government of the Hong Kong SAR. 
FSNL and JPM acknowledge financial support of the 
Fundagao para a Ciencia e Tecnologia through the grants 
PTDC/FIS/102742/2008 and CERN/FP/109381/2009. 



Acknowledgments 



NMG acknowledges a postdoctoral fellowship from 
CONACYT-Mexico. The work of TH is supported by an 



[1] S. Perlmutter et ai, Astrophys. J. 517, 565 (1999); 
A. G. Riess et al, Astron. J. 116, 1009 (1998); 
A. G. Riess et ai, Astrophys. J. 607, 665 (2004); A. 
Grant et al, Astrophys. J. 560 49-71 (2001); S. Perlmut- 
ter, M. S. Turner and M. White, Phys. Rev. Lett. 83 670- 
673 (1999); C. L. Bennett et al, Astrophys. J. Suppl. 148 
1 (2003); G. Hinshaw et al, | arXiv:astro-ph/0302217] ; 

E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. 
Phys. D 15, 1753 (2006). 

[2] K. Koyama, Gen. Rel. Gray . 40, 421 (2008); 

F. S. N. Lobo, larXiv:0807.T640l [gr-qc]; A. De Felice 
and S. Tsujikawa, Living Rev. Rel. 13, 3 (2010). H. A. 
Buchdahl, Mon. Not. Roy. Astron. Soc. 150, 1 (1970); 
R. Kerner, Gen. Rel. Grav. 14, 453 (1982); J. P. Du- 
ruisseau, R. Kerner and P. Eysseric, Gen. Rel. Grav. 
15, 797 (1983); J. D. Barrow and A. C. Ottewill, J. 
Phys. A: Math. Gen. 16, 2757 (1983); A. A. Starobinsky, 
Phys. Lett. B 91, 99 (1980); S. M. Carroll, V. Duvvuri, 
M. Trodden and M. S. Turner, "Is cosmic speed-up due 
to new gravitational physics?," Phys. Rev. D 70, 043528 

(2004) ; L. Amendola, D. Polarski and S. Tsujikawa, Phys. 
Rev. Lett. 98, 131302 (2007); S. Capozziello, S. Nojiri, 
S. D. Odintsov and A. Troisi, Phys. Lett. B639, 135 

(2006) ; S. Nojiri and S. D. Odintsov, Phys. Rev. D74, 
086005 (2006); M. Amarzguioui, O. Elgaroy, D. F. Mota 
and T. Multamaki, Astron. Astrophys. 454, 707 (2006); 
L. Amendola, R. Gannouji, D. Polarski and S. Tsu- 
jikawa, Phys. Rev. D75, 083504 (2007); T. Koivisto, 
Phys. Rev. D 76, 043527 (2007); A. A. Starobinsky, 
JETP Lett. 86, 157 (2007); B. Li, J. D. Barrow and 

D. F. Mota, Phys. Rev. D 76, 044027 (2007); S. E. Perez 
Bergliaffa, Phys. Lett. B642, 311 (2006); G. Cognola, 

E. Elizalde, S. Nojiri, S. D. Odintsov and S. Zerbini, 
JCAP 0502, 010 (2005); V. Faraoni, Phys. Rev. D72, 
061501 (2005); V. Faraoni, Phys. Rev. D72, 124005 

(2005) ; L. M. Sokolowski, Class. Quant. Grav. 24, 3391 

(2007) ; G. Cognola, M. Gastaldi and S. Zerbini, Int. J. 
Theor. Phys. 47, 898 (2008); C. G. Bohmer, L. Hol- 
lenstein and F. S. N. Lobo, Phys. Rev. D 76, 084005 
(2007); S. Carloni, P. K. S. Dunsby and A. Troisi, Phys. 
Rev. D 77, 024024 (2008); K. N. Ananda, S. Carloni 
and P. K. S. Dunsby, Phys. Rev. D 77, 024033 (2008); 
S. Capozziello, R. Cianci, C. Stornaiolo and S. Vignolo, 
Class. Quant. Grav. 24, 6417 (2007); S. Tsujikawa, Phys. 
Rev. D 77, 023507 (2008); S. Nojiri, S. D. Odintsov 
and P. V. Tretyakov, Phys. Lett. B 651, 224 (2007); 



S. Nojiri and S. D. Odintsov, Phys. Lett. B 652, 343 

(2007) ; G. Cognola, E. Elizalde, S. Nojiri, S. D. Odintsov, 
L. Sebastiani and S. Zerbini, Phys. Rev. D 77, 046009 

(2008) ; T. Chiba, Phys. Lett. B575, 1 (2003); A. L. Er- 
ickcek, T. L. Smith and M. Kamionkowski, Phys. Rev. 
D74, 121501 (2006); T. Chiba, T. L. Smith and 
A. L. Erickcek, Phys. Rev. D75, 124014 (2007). S. No- 
jiri and S. D. Odintsov, Phys. Lett. B 659, 821 (2008); 
S. Capozziello, A. Stabile and A. Troisi, Phys. Rev. D 76, 
104019 (2007); S. Capozziello, A. Stabile and A. Troisi, 
Class. Quant. Grav. 25, 085004 (2008); G. J. Olmo, Phys. 
Rev. D75, 023511 (2007); W. Hu and I. Sawicki, Phys. 
Rev. D76, 064004 (2007); S. Nojiri and S. D. Odintsov, 
Phys. Rev. D68, 123512 (2003); V. Faraoni, Phys. 
Rev. D74, 023529 (2006); T. Faulkner, M. Tegmark, 

E. F. Bunn and Y. Mao, Phys. Rev. D76, 063505 

(2007) ; P. J. Zhang, Phys. Rev. D 76, 024007 (2007); 
S. Capozziello and S. Tsujikawa, Phys. Rev. D 77, 
107501 (2008); I. Sawicki and W. Hu, Phys. Rev. D75, 
127502 (2007); L. Amendola and S. Tsujikawa, Phys. 
Lett. B 660, 125 (2008); S. Capozziello, V. F. Cardone 
and A. Troisi, JCAP 0608, 001 (2006); S. Capozziello, 
V. F. Cardone and A. Troisi, Mon. Not. R. Astron. 
Soc. 375, 1423 (2007); A. Borowiec, W. Godlowski and 
M. Szydlowski, Int. J. Geom. Meth. Mod. Phys. 4 (2007) 
183; C. F. Martins and P. Salucci, Mon. Not. Roy. As- 
tron. Soc. 381, 1103 (2007); C. G. Boehmer, T. Harko 
and F. S. N. Lobo, Astropart. Phys. 29, 386 (2008); 
C. G. Boehmer, T. Harko and F. S. N. Lobo, JCAP 0803, 
024 (2008); O. Bertolami, C. G. Boehmer, T. Harko 
and F. S. N. Lobo, Phys. Rev. D 75, 104016 (2007); 
O. Bertolami and J. Paramos, Phys. Rev. D 77, 084018 

(2008) ; O. Bertolami, F. S. N. Lobo and J. Paramos, 
Phys. Rev. D 78, 064036 (2008); arXiv:0811.2876l [gr- 
qc]. F. S. N. Lobo and M. A. Oliveira, Phys. Rev. 
D 80, 104012 (2009); C. G. Boehmer, L. Hollenstein, 

F. S. N. Lobo and S. S. Seahra. larXiv:1001. 12661 [gr-qc]: 
N. M. Garcia and F. S. N. Lobo, Phys. Rev. D 82, 104018 
(2010); N. M. Garcia and F. S. N. Lobo, Class. Quant. 
Grav. 28, 085018 (2011); T. Harko and F. S. N. Lobo, 
larXiv: 1007.44151 [gr-qc]; T. Harko and F. S. N. Lobo, 
Eur. Phys. J. C 70, 373 (2010); T. Har ko, F. S N. Lobo, 
S. Nojiri and S. D. Odintsov. larXiv:1104 .2669 [gr-qc]. 

[3] H. Weyl, Space, Time, Matter, Chapter IV, New York, 
Dover (1921); A. Eddington, The Mathematical Theory 
of Relativity, Chapter IV, London, CUP (1924); K. Lane- 



9 



zos, Ann. Math. 39, 842 (1938); H. A. Buchdahl, Proc. 

Edin. Math. Soc. 8, 89 (1948); A. Pais and G. E. Uh- 

lenbeck, Phys. Rev. 79, 145 (1950); R. Utiyama and B. 

de Witt, J. Math. Phys. 3, 608 (1962); P. Havas, Gen. 

Rel. Grav. 8, 631 (1977); K. Stelle, Gen. Rel. Grav. 5, 

353 (1978); P. D. Mannheim and D. Kazanas, Astroph. 

Journ. 342, 635 (1989); D. Kazanas and P. D. Mannheim, 

Astroph. Journ. Supp. Series 76, 421 (1991). 
[4] Nojiri, S., Odintsov, S.D., and Sasaki, M. 2005, Phys. 

Rev. D, 71, 123509. 
[5] Nojiri, S., and Odintsov, S.D 2007. Int. J. Geom. Meth. 

Mod. Phys., 4, 115; S. Nojiri and S. D. Odintsov, 

larXiv:1011.0544l [gr-qc] . 
[6] Cognola, G., Elizalde, E., Nojiri, S., Odintsov, S.D., and 

Zerbini, S. 2006, Phys. Rev. D, 73, 084007; Nojiri, S., 

Odintsov, S.D., and Gorbunova, O.G. 2006, J. Phys. A: 

Math. Gen., 39, 6627. 
[7] Sotiriou, T.P. 2007, arXiv:gr-qc/0710.4438. 
[8] S. Nojiri and S. D. Odintsov, Phys. Lett. B 631, 1 (2005). 
[9] S. Nojiri, S. D. Odintsov and M. Sami, Phys. Rev. D 

74, 046004 (2006); S. Nojiri and S. D. Odintsov, J. Phys. 

Conf. Ser. 66, 012005 (2007); G. Cognola, E. Elizalde, 

S. Nojiri, S. Odintsov and S. Zerbini, Phys. Rev. D 75, 

086002 (2007); B. Li, J. D. Barrow and D. F. Mota, Phys. 

Rev. D 76, 044027 (2007); M. Gurses, Gen. Rel. Grav. 

40, 1825 (2008); A. De Felice and S. Tsujikawa, Phys. 
Lett. B 675, 1 (2009); M. Alimohammadi and A. Ghalee, 
Phys. Rev. D 79, 063006 (2009); C. G. Boehmer and 
F. S. N. Lobo, Phys. Rev. D 79, 067504 (2009); K. Ud- 
din, J. E. Lidsey and R. Tavakol, Gen. Rel. Grav. 

41, 2725 (2009); S. Y. Zhou, E. J. Copeland and 
P. M. Saffin, JCAP 0907, 009 (2009); A. De Felice 
and T. Suyama, JCAP 0906, 034 (2009); N. Goheer, 
R. Goswami, P. K. S. Dunsby and K. Ananda, Phys. 
Rev. D 79, 121301 (2009); J. Sadeghi, M. R. Setare and 
A. Banijamali, Phys. Lett. B 679, 302 (2009); J. Sadeghi, 
M. R. Setare and A. Banijamali, Eur. Phys. J. C 64, 433 



(2009); A. De Felice and S. Tsujikawa, Phys. Rev. D 80, 
063516 (2009); M. Alimohammadi and A. Ghalee, Phys. 
Rev. D 80, 043006 (2009); K. Bamba, C. Q. Ceng, S. No- 
jiri and S. D. Odintsov, Europhys. Lett. 89, 50003 (2010); 
M. Mohseni, Phys. Lett. B 682, 89 (2009); D. Bazeia, 
B. Carneiro da Cunha, R. Menezes and A. Y. Petrov, 
Phys. Lett. B 649, 445 (2007) D. Bazeia, R. Menezes 
and A. Y. Petrov, Eur. Phys. J. C 58, 171 (2008). 

[10] G. Cognola, M. Gastaldi and S. Zerbini, Int. J. Theor. 
Phys. 47, 898 (2008). 

[11] S. Nojiri, S. D. Odintsov and P. V. Tretyakov, Prog. 
Theor. Phys. Suppl. 172, 81 (2008). 

[12] K. Bamba, S. D. Odintsov, L. Sebastiani and S. Zerbini, 
Eur. Phys. J. C 67, 295 (2010). 

[13] S. W. Hawking and G.F.R. Ellis, The Large Scale Struc- 
ture of Spacetvme, (Cambridge University Press, England, 
1973). 

[14] J.H. Kung, Phys. Rev. D 52, (1995) 6922; Phys. Rev. 
D 53, (1996) 3017; S.E. Perez Bergliaffa, Phys. Lett. B 
642, (2006) 311; S. Carroll, Spacetime and Geometry: 
An Introduction to General Relativity, (Addison Wesley, 
New York, 2004); J. Santos and J.S. Alcaniz, Phys. Lett. 
B 619, (2005) 11; M. Visser, Science 276, (1997) 88; 
Phys. Rev. D 56, (1997) 7578; J. Santos, J.S. Alcaniz 
and M.J. Reboucas, Phys. Rev. D 74, (2006) 067301; J. 
Santos, J.S. Alcaniz, N. Pires and M.J. Reboucas, Phys. 
Rev. D 75, (2007) 083523; A. A. Sen and R. J. Scherrer, 
Phys. Lett. B 659, 457 (2008); J. Santos, J. S. Alcaniz, 
M. J. Reboucas and N. Pires, Phys. Rev. D 76, 043519 
(2007); Y. Gong and A. Wang, Phys. Lett. B 652, 63 
(2007). 

[15] J. Santos, J. S. Alcaniz, M. J. Reboucas and F. C. Car- 

valho, Phys. Rev. D 76, 083513 (2007). 
[16] D. Rapetti, S.W. Allen, M.A. Amin and R.D. Blandford, 

Mont. Not. R. Soc. 375, (2007) 1510. 
[17] N. J. Poplawski, Class. Quant. Grav. 24, 3013 (2007). 



